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ABSTRACT 

We present the results of Hi 21-cm line observations to explore the nature of 
the high- velocity (HV) Hi gas at i ~ 173°. In low- resolution Hi surveys this 
HV gas appears as faint, wing-like, Hi emission that extends to velocities beyond 
those allowed by Galactic rotation. We designate this feature FVW (Forbidden 
Velocity Wing) 172.8+1.5. Our high-resolution (3/4) Arecibo Hi observations 
show that FVW 172.8+1.5 is composed of knots, filaments, and ring-like struc- 
tures distributed over an area a few degrees in extent. These HV Hi emission 
features are confined within the limits of the HII complex G173+1.5, which is 
composed of five Sharpless HII regions distributed along a radio continuum loop 
of size 4°. 4 x 3°. 4, or ~ 138 pc x 107 pc, at a distance of 1.8 kpc. G173+1.5 
is one of the largest star-forming regions in the outer Galaxy. We demonstrate 
that the HV Hi gas is well correlated with the radio continuum loop and that 
the two seem to trace an expanding shell. The expansion velocity of the shell 
is large (55 km s _1 ) suggesting that it represents a supernova-remnant (SNR). 
We derive physical parameters for the shell and show these to be consistent with 
the object being a SNR. We also detect hot X-ray emitting gas inside the HII 
complex by analyzing the ROSAT all-sky X-ray background survey data. This 
also supports the SNR interpretation. We conclude that the HV Hi gas and the 
X-rays are most likely the products of a supernova explosion(s) within the HII 
complex, possibly in a cluster that triggered the formation of these HII regions. 



Subject headings: supernova remnants — ISM: individual (FVW172.8+1.5) - 
radio lines: ISM — ISM: Hll regions - stars: formation 
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Introduction 



In large-scale (£,v) diagrams of Hi 21-cm line emission in the Galactic plane, there are 
many small, faint, high- velocity "wing-like" features exte nding to velocities w ell beyond the 
maximum or minimum permitted by Galactic rotation fjKoo fc Kangi l2004j ) . These "For- 
bidden Velocity Wings (FVWs)" are likely due to energetic phenomena in the Galaxy, as 
they are confined to small areas (5, 2°), and projec t smoothly beyond the general Galac- 
tic emission to high velocities. iKoo fc Kangi (120041 ) suggested that some of these FVWs 
may represent the expanding shells of "missing" supernova remnants (SNRs), which are 
not included in existing SNR catalogs. The basic idea is that old SNRs are faint in the 
radio continuum, making it difficult to identify them because of both the confusion due to 
the Galactic b ackgr ound emission, and observational limitations (e.g., Brogan et al. 2006). 
Koo fc Kangi (120041 ) considered the fact that an old SNR should possess a fast-expanding 
Hi shell that will still be present as a coherent entity after the remnant becomes too faint 
to be visible in the radio continuum. If its expansion velocity is greater than the minimum 
or maximum velocities permitted by Galactic rotation, then an old SNR shell, o r part of it 



(mayb e the caps), could be detected as high- velocity Hi gas, e.g., FVWs. Indeed, IKoo et al. 



( 120061 ) have carried out high-resolution Hi line observations toward one FVW and detected 
a rapidly expanding (80 km s _1 ) Hi shell, the parameters of which are consistent with those 
of the remnant of a SN that exploded some 0.3 Myr ago. 



Rece ntly. iKang fc Kool (2007) id entified 87 FVWs from the Leiden/Dwingeloo Hi sur- 
vey (LPS; lHartmann fe Burtonlll997l ) and the Southern Galactic Plane Survey data (SGPS; 
McClure- Griffiths et al.ll2005l ). Among these, 6 FVWs are found to be coincident with SNRs, 
4 with nearby galaxies, and 3 with high-velocity clouds. The rest (85%) are not associated 
with any obvious objects that could be responsible for their high velocities. We have since 
been making follow-up high-resolution Hi observations of these FVWs of unknown nature 
using the Arecibo 305-m and Green Bank 100-m radio telescopes in order to identify their 
natures and origins, and have now observed ~ 30% of them. We find that about 40% of 
the observed FVWs have shell-like morphologies. More than half of these are apparently 
expanding at ^ 50 km s _1 , which supports the possibility of some FVWs being old SNRs. 
The other 60% show irregular structures consisting of filaments and clumps. Most of these 
show faint bumps in their line profiles indicating high-velocity Hi clouds. They might cor- 
respond to the fast-moving, compact clouds in the disk, or in the disk-halo interface, that 



have been recently detected in sensitive, high-resolution Hi surveys (jStanimirovic et al.ll2006 



Begum et al.ll2010l and references therein). The full results from our survey will be presented 
in a separate paper. 



Here, we report Arecibo Hi 21-cm line observations of FVW172.8+1.5. This particular 
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FVW is unique in that it is associated with an HII complex in the outer Galaxy. We show 
that FVW 172.8+1 .5 is most likely an old SNR produced in this complex. The catalog of 



Kang fc Kool (120071 ) lists this feature as two objects, FVW173.0+0.0 and FVW173.0+3.0. 



However, the Arecibo Hi image shows that these are likely to be parts of a single coherent 
object, which we will call FVW172.8+1.5. In Section 2, we describe the Hi observations. 
The Hi results are presented in Section 3. In Section 4, we present a multi-wavelength view 
of the HII complex in this region and investigate its association with FVW 172.8+1.5. We 
discuss the origin of FVW172.8+1.5 and the star-formation history in this area in Section 
5, presenting a summary in Section 6. 



2. The Hi 21-cm Line Observations 



Hi 21-cm line observations of FVW172.8+1.5 were made in 2006 October with the 
Arecibo 305-m telescope using the 7-beam, dual-polarization, Arecibo L-band Feed Array 
(ALFA). The beam-to-beam spacing for the scanning pattern used was 1'.8. GALSPECT, a 
dedicated spectrometer for Galactic Hi ALF A surveys, was used as the backend. Both wide- 
and narrow-band spectra were acquired (see IStanimirovic et aLj|2006j , for details). Here we 
present data from the narrower 7.14-MHz band, which is analyzed into 8192 equally spaced 
channels, giving an unsmoothed velocity resolution of 0.18 km s _1 . 

The field observed covers an area of 7?5 x4?5 centered at (a J) = (5 h 37 m , 35°40') (J2000). 



The observations employed the basket- weaving technique (IPeek fc Heilesl l2008t iPeek et al. 



20071 ). which scans the sky with a zig-zag pattern by driving the telescope up and down on the 
prime meridian at a speed of 1.5 arcmin /sec. For bandpass correction, a "Least-Squares Fre- 
quency Switching" (LSFS; iHeilesI 120071 ) calibration was perform ed each day. Data re duction 
used the IDL pipe-line codes developed by the Berkeley group (IPeek fc Heilesl 120081 ). Using 
these, the data are converted into a brightness temperature cube. The cont amination due 
to the coma sidelobes of the ALFA off-centre beams has been corrected (see Peek & Heiles 



20081 . for the details of data reduction). The final cube has a velocity resolution of 0.74 km 



a spatial HPBW of 3.4 arcmin, and an rms of 0.13 K. 



3. The Hi Results 

Fig. [1] shows channel images of the Arecibo Hi data. The Hi emission associated with 
FVW172.8+1.5 is visible for +19.5 km s _1 < v LS r < +50.4 km s -1 . At the highest velocities 
(+40.1 km s -1 < t>LSR ~ +50.4 km s" 1 ), the Hi appears to be largely separated into two con- 
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centrations centered at ~ (5 h 37 m , +36°30') and ~ (5 h 27 m , +34°30'), These were previously 
designated FVW173. 0+3.0 and FVW173. 0+0.0 respectively. The northern concentration 
(FVW173. 0+3.0) looks diffuse and clumpy, and appears to form a ring-like structure with 
a diameter of ~ 2°. In contrast, the southern concentration (FVW173. 0+0.0) appears to 
be confined to a few small (< 30') regions. At lower velocities (+16.6 km s -1 < ^lsr < 
+35.0 km s _1 ), a thin, filamentary feature that surrounds the two concentrations appears, 
so that the entire structure appears to be a large shell with an extent of 4°. 4 x 3.4°. The 
crosses in Fig. [T] and the ellipse superposed on the +24.7 km s _1 channel map mark the 
approximate center and boundary of the shell structure. At +14.4 km s" 1 and lower ve- 
locities, the Hi emission from the Galactic ISM dominates, so that no particular Hi feature 
associated with the source is recognizable (see also § 4.1.1). 

The various Hi features described above are also visible in Fig. [2J which is a three-color 
image of the Hi gas integrated over different velocit y intervals. In th is figure, we super- 
pose the Effelsberg 11-cm radio continuum contours (jFiirst et al.lll990l ). It is seen that the 
high-velocity Hi emission features lie essentially within the radio continuum filaments. The 
northwestern high velocity Hi gas, represented by red and green colors, are confined within 
the northern part of the continuum structure, while the Hi features at lower velocities 
(colored blue) lie along the inner boundary of the outer radio continuum filaments. The 
morphological relation between the Hi and the radio filaments strongly suggests their as- 
sociation. In addition, there is little possibility of chance alignment of the two structures, 
because confusion is low along the line of sight in the direction of I ~ 173°. The increase of 
the surface area and/or the size of the Hi -emitting region with decreasing velocity suggests 
that we are looking at the receding portion of an expanding shell. The Hi morphology does 
not exactly match that of a uniformly expanding spherical shell, but this could be due to a 
non-uniform and inhomogeneous ambient interstellar medium. The radio continuum struc- 
ture is part of the HII complex in this area, and we will investigate the association of the 
two at other wavelengths in Section 4 and their physical properties and origin in Section 5. 

The velocity structure of the high- velocity Hi features is shown in Fig. [3], which shows the 
peak Hi brightness temperature in Galactic longitude as a function of Galactic latitude and 
velocity. This position-velocity diagram shows the high- velocity gas over the entire field in a 
single diagram and enhances the appearance of small-scale wing-like features. Fig. [3] shows 
that FVW172.8+1.5 is composed of many small-scale, high-velocity wings spatially confined 
to — 1° < b < +4°, which suggests an inhomogeneus nature for the shell. These wings are 
mostly straight but some have the shape of an extended ring, e.g., that at b = —0.2°. This 
indicates that the shell is composed of small clumps. Since only the high- velocity end portions 
of the wings are visible, their central velocities are usually not accessible. However, a few line 
profiles in the area near (5 h 34.5 m , +36°03') show discrete spectral features at high velocities 
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(Fig. H]). A Gaussian fit to these features resuits in a central velocity of +33 — [-35 km s~ x 
and a velocity dispersion of 4.3 - 6.4 km s _1 , corresponding to kinetic temperatures of 2,200 
- 4,900 K if the broadening is entirely due to thermal motions. 

4. The i = 173° region in Multi Wavebands 
4.1. Radio Continuum 



4.1.1. The Ell Complex Gl 73+1.5 



The region of Auriga within which FVW172.8+1.5 lies is an HII complex composed of 
Sharpless HII regions and OB associations in the Perseus arm (Fig. [5]). The H II regions are 
organ ized along a large (~ 7° x 4°) filamentary structure r esembling a bow tie (jrCerton et al. 



20071 ). The Hll regions have been studied previously (e.g. Ilsrael &: Fellilll978l ). and they are 
largely at two different distances (see Tabled]). To the northeast of the structure, S231, S232, 
S233, a nd S235 are known to be associated with a giant molecular cloud at a distance of 



1.8 kp c (lEvans fc BlaiJl981i l. with CO radial velocities of -18.1 to -23.0 km s" 1 (Blitz et al. 



19821 ). in which active star formation is on-going (see § 4.3). The distances determined by 
spectrophotometry of the exciting stars of individual HII regions range from 1.0 to 2.3 kpc, 
but the distance to the associated molecular cloud (1.8 kpc) is often adopted as the distance 
to the whole HII complex. To the south of the complex, two l°-sized HII regions S229 (IC 
405) and S236 (IC 410) and the large (~ 5°) diffuse HII region S230 are present. S236 is 
associated with a molecular cloud at —7.2 km s -1 and might be at a considerably greater 
distance bec ause the estimated distance to its central star cluster (NGC 1893) ranges from 
3.2 to 6 kpc (jSharma et al.ll2007l . and references th erein). For the ot her two HII regions (S229 
and S230), not much has been known, although iFich et al.l (jl990l ) measured Ha velocities 
of +4.4 and 0.0 km s^ 1 toward S229 and S230, respectively. Toward S229, CO J = 1-0 
emission at +6.7 km s -1 has been detected, but it may not be related to the HII region 



(IBlitz et al.lll982l ). The small HII regions S234 (IC 417) and S237 in the middle of the "bow- 
tie" have velocities of —13 and —4 km s _1 , similar to the first and second groups respectively. 
A recent photometric determination of the distance to the central cluster (Stock 8) of S234 
also agrees with the distance to the northwestern group of HII regions, i.e., 2.05 ± 0.10 kpc 
( IJose et al.ll2008l ). We therefore consider that the Sharpless HII regions S231-S235 form the 
active star formation complex G173+1.5 at a distance of 1.8 kpc at t>LSR ~ ~~ 20 km s _1 . 
The other HII regions have considerably different velocities and distances, and they are not 
considered to be associated with this complex. 



The morphology of the region suggests an association between the large radio continuum 
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filamentary structures and the HII regions. Filaments A through D (see Fig. 5) have similar 
concave shapes suggesting that they have resulted from a common source within the HII 
complex. For the radio filament A, we have found an HI filamentary structure along the 
radio feature at t>LSR =—25 to —28 km s _1 (Fig. E]). The correlation is not perfect in the 
sense that the HI feature looks less curved than the continuum filament and seems to extend 
outside of the field shown in Fig. 6. However, the velocity is close to the CO velocity of 
S232, and the two could be associated. We could not find any Hi structures associated with 
the other continuum filaments. Filament E has an enhanced brightness and its curvature is 
opposite to those of Filaments A-D. It is likely that the enhanced brightness and the convex 
shape is due to interaction with a dense ambient medi um there, althoug h no responsible 
molecular cloud is seen in the available CO survey data ( jDame et al.ll200ll see also Fig. 9). 
Hence, filaments A-E are likely to be associated with the HII complex G173+1.5, although 
we have only circumstantial evidence, except perhaps for filament A. Kinematic evidence is 
needed to confirm the association. 



4- 1.2. The Radio Continuum Spectrum of Filament A 
The nature of the radio emission from the l arge continuum filam entary structure has 



been discussed previously, although not in detail. iKerton et al.l ( 120071 ) noted this structure 
in the 1420 MHz Canadian Galactic Plane Survey map, and argued that the filaments D & E 
are thermal emission based on either the presence of corresponding i nfrared emission and/or 



the rising spectrum between 408 and 1420 MHz. iGao et al.l (120101 ) claim that filament A 



shows a thermal spectrum between the higher frequencies of 1.4 and 5 GHz. 

Using radio continuum data from 325 MHz to 2.7 GHz available on the web (Table. 
we have attempted a study of the continuum spectrum of the filaments associated with the 
structure. The faint filamentary features are visible at most of these frequencies. Here we 
will focus on Filament A at (£, b) ~ (172?5, +3?5), which stands out at 325 MHz in WENSS. 
To estimate values of spectral index, we (1) converted the WENSS data to brightness tem- 
perature, thus bringing all surveys to the same intensity units, (2) convolved all images to 
have a circular beam of FWHM = 4'. 7, the declination beam size of the CGPS 408 MHz 
image at 5 = 37°, (3) gridded the images onto the 2'-spaced pixels of the Effelsberg data, and 
(4) subtracted out all point sources. The final images are shown in Fig. [7J where the dotted 
and solid boxes mark Filament A and the Hll region, S232, which provides a reference with 
a thermal spectrum. Finally, we removed a smooth, linear, large-scale background emission 
gradient using two neighboring areas, marked by pairs of boxes flanking the Hll region and 
the filament in Fig. [71 
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The continuum spectra of Filament A and S232 are shown in Fig. [8] (left) by filled and 
open circles respectively. The total flux densities are given in Table [31 The error bars in Fig. [8] 
represent standard deviations in the neighboring areas after removal of smooth backgrounds. 
Using all flux density values from 325 to 2695 MHz, the estimated spectral index (S u oc v~ a ) 
of S232 is a = -0.09 ± 0.02, while that of Filament A is a = +0.23 ± 0.06. We also derived 



spectral indices for Filament A and S232 using T-T plots ( [Turtle et al.lll962l ). The regions 
used for making the T-T plots are those over which the total flux densities were derived. 
Fig. El (right) shows a T-T plot between 325 and 2695 MHz of Filament A. The points are 
well described by a linear fit. The slopes of the plots, and the values of a derived between 
different pairs of frequencies, are summarized in Table HI For S232, the derived spectral 
indices are consistent, ranging from a = —0.04 to +0.09 which, together with the total flux 
density spectral index of a = —0.09, is reasonably close to the value of a = +0.1 expected 
for an optically thin thermal source. However, the spectral indices derived for Filament A 
are not straightforward to interpret. The frequency pairs that include 325 MHz result in 
a = +0.50 ± 0.11 and +0.32 ± 0.06, consistent with non-thermal emission, while the pairs 
that include 408 MHz yield a = +0.07 ± 0.10 and -0.04 ± 0.08, compatible with thermal 
emission. These trends are consistent with the plot of the total flux densities (Fig. [5] left). 
In that plot, if either the 325 MHz or the 408 MHz data point is excluded, then the spectral 
index of Filament A would respectively be flatter or steeper than the a = +0.23 ± 0.06 
derived using data at all frequencies. 

As mentioned at the beginning of this section, previous studies claimed that the radio 
emission from some of these filaments is of thermal origin. However, our analysis in this 
section indicates the presence of a non-thermal component in filament A. The fact that the 
spectral index of S232, as derived from the same data, is consistent with optically-thin ther- 
mal emission supports our result. Therefore, it is more likely that the filament has both 
thermal and non-thermal components (see also § 4.3). If thermal and non-thermal compo- 
nents co-exist in the filament, a steep spectrum at lower frequencies, and a flat spectrum 
at higher frequencies (as found by Gao et al.), might be expected. Also, if thermal and 
non-thermal components co-exist in the continuum emission of Filament A, this could be 
the case for Filaments B-E. We have examined the polarization maps of this area. The G 173 



comp l ex region was co vered by the polarization surveys at 1.4 and 5 GHz (ILandecker et al. 



20101 ; iGao et al.l 120101 ) . While the 1.4 GHz polarization data of Landecker et al. do not 
show any apparent counterparts to the G173 structure, the 5 GHz polarization data of Gao 
et al. do show some filamentary features associated with the structure. In the polarized 
intensity map of Gao et al. (their Fig. 17), two depolarized filaments, indicating a thermal 
nature, are recognizable; one at (£,b) ~ (174.3, +2.0), which corresponds to Filament B, 
and the other at ~ (171.0, +2.5), located between Filaments A and D. Other than B, the 
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continuum filaments have no apparent counterparts in this image. Hence, it is difficult to 
discuss the nature of Filaments A - E using the current polarization data, with the exception 
that a thermal component appears to dominate in Filament B. The nature of the continuum 
filaments needs to be investigated with more sensitive full-Stokes continuum data. 



4.2. X-rays 



In the ROSAT Survey Diffuse X-ray Background Map (ISnowden et al.lll997l ). there is 
faint, extended X-ray emission projected against the HII complex G173+1.5. (Fig. O top 
frames). The emission is mainly from two regions, one within the S231-235 complex and 
the other associated with S229/S236. The emission is not visible in the soft band (1/4 keV) 
but appears at the hard bands (3/4 and 1.5 keV). This is consistent with the emission being 
associated with G173+1.5 because if we a dopt Ay = 1.2 m ag (or N{H) = 2.3 x 10 21 cm -2 ) 



corresponding to the extinction to S234 (IJose et al.l 120081 ) as the extinction to G173+1.5, 
the transmission of 1/4 keV ph otons is essen tially zero ( ^ 10%) while it is 40% and 80% at 
3/4 and 1.5 keV, respectively ( jSewardl l2000f ) . The emission is not uniform and the ratio of 
1.5 to 3/4 keV intensities vary over the field. Note that, since the contributions from point 
sources have been removed from the map, the emission is from diffuse sources, although 
some compact sources might have been included, e.g., the bright spot near S234. For the 
present analysis, we simply derive the X-ray photon counts inside the solid ellipse marked in 
Figj9j The background is estimated from the surrounding regions marked by the two dotted 
ellipses. The derived count rates over 8.8 deg 2 are 0.76 ±0.08 cts s _1 and 0.78 ±0.11 cts s _1 
in 3/4 keV and 1. 5 keV bands, respecti vely. The ratio of 1.5 keV to 3/4 keV is 1.0 ± 0.2. 
Using the model of ISnowden et al.l (119971 ). the ratio corresponds to either a power-law photon 
spectrum (E~ a , where E is energy) of index a ~ 4 or thermal spectrum with temperature 
of ~ 10 6 ' 9 K for N(H) = 2.3 x 10 21 cm~ 2 . The power-law index 4 is considerably steeper 
than that of a pulsar wind nebula, e.g., ~ 2 for the Crab nebula, so that we may rule 
out a non-the rmal origin. Th e diffu se extended morphology also supports a thermal origin. 
According to ISnowden et al.l ( 119971 ) , the attenuated count rate of thermal X-ray emission 
from 10 6 ' 9 K gas at 3/4 keV is ~ 0.01 cts s _1 acrmin -2 for an emission measure of 1 cm -6 
pc. Therefore, the observed average 3/4 keV-band count rate of 2.4 x 10 -5 cts s _1 arcmin -2 
implies an average emission measure of 0.0024 cm -6 pc. If we use 4R/3, where R = 60 pc 
is the geometrical mean radius of the shell, as the mean depth along the line of sight, the 
mean electron density is about 5.5 x 10~ 3 cm~ 3 . This implies a thermal energy for the X-ray 
emitting gas of ~ 3 x 10 50 ergs, assuming that it fills the entire interior of the shell. If instead 
the X-ray emitting gas fills only parts of the shell interior, then the thermal energy will be 
somewhat smaller than this. The derived thermal energy is close to that of an old SNR, 
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which, together with the high temperature, suggests that the hot gas might have originated 
from an SN explosion within the HII complex G173+1.5. 



4.3. Ha and CO 



The radio continuum structure of G173+1.5 has counterparts in Ha . Fig. |H] (bottom left) 
is the Ha full-sky map (6' FWHM resolution) composite of the Virginia Tech Spectral line 
Survey (VTSS), the Southe rn Ha Sky Survey Atlas (SHASSA), and the Wisconsin H-Alpha 
Mapper ( WHAM) survey (IFinkbeinerl 120031 ). In the northern area covered by the VTSS 



Ha image (foennison et all boooh delicate Ha filaments are clearly visible, well correlated 



with the radio continuum filaments. If we compare in detail, however, there is a slight shift 
between the radio and Ha filaments; the Ha emission peaks at slightly lower latitudes, i.e. 
slightly inside, the complex (Fig. [TO]) . If the radio feature is due to thermal emission, the 
ratio of the two intensities are roug h ly co nstant as both depend on the emission measure. 
Using the formulae given in ISpitzerl (119781 ). the Ha intensity from ionized gas at 10,000 K 



is given by I (Ha) ~ 0.36(n p /n e )EM Rayleigh (R) where EM = J n 2 e ds (cm -6 pc), and n p 
and n e are proton and electron densities, resp ectively . The corresponding radio brightness 
temperature at 11 cm, using the formulae from lDraind (120111 ). is T& ~ 0.40 x 10~ 3 (n p /n e )EM 
K. Therefore, we have T^/ I(Ha) ~ 4 x 10~ 3 K R _1 adopting an extinction of Ay = 1.4 mag. 
At T = 5, 000 K, this ratio would be slightly lower. Since the Ha intensity is normalized by 
20 R and the 11-cm brightness is normalized by 0.1 K, the thermal 11-cm emission would 
have a comparable (~ 0.8) normalized intensity in Fig. [10J In fact, Fig. [10] shows that at 
the Ha peak positions this is indeed the case. On the other hand, at the radio continuum 
maxima, the 11-cm brightness is considerably greater than the Ha intensity, which implies 
a non-thermal origin. Therefore, the 11-cm continuum emission appears to be composed of 
both thermal and non-thermal components, which is consistent with our conclusion from the 
previous section. 

Figure [§ ] (bottom right) sh ow the CO J=l-0 intensity map integrated over t> L sR = —25 to 
—15 km s _1 (jDame et al.ll200ll ). There is one giant molecular cloud at (£, b) ~ (173. °5, 2.°5). 
Four Hll regions, S231, S232, S233, and 235, are clustered around this molecular cloud. 
CO observations of this cloud s how close position al correlations between the optical Hll 
regions and the molecular cloud (jHeyer et al.lll996l ). The giant molecular cloud appears as 
interconnected filaments with bright CO emission located at the projected optical edges of 
S235 and S232. S233 lies within a void in the CO emission, implying complete ionization and 
photodissociation of the molecular material. It has been known that active star formation 
is ongoing around these Hll regions, and that jets and outflows from young stellar objects 



have been observed in molecular lines (e.g., IShepherd fc Watson 



2002 



Beuther et al 



2002|). 



5. Discussion 

5.1. The Origin of FVW172.8+1.5 

Our results indicate that the fast-moving Hi gas in FVW172.8+1.5 is confined within the 
HII complex G173+1.5 and has good spatial correlation with the radio continuum structure. 
This strongly suggests its origin lies within the complex. Morphological association between 
the Hi emission and the continuum/Ha filaments further suggests that they probably trace 
the same object, i.e., an expanding shell. The X-rays detected within the complex may 
also be associated with this shell. In this section, we first derive the parameters that an 
expanding shell would possess, and then explore its possible origin. We take the systemic 
velocity and distance of both the HII complex and FVW172.8+1.5 to be t> sys ~ —20 km s _1 
and 1.8 kpc. 



5.1.1. The parameters of FVW172.8+1.5 as an expanding shell 



The total extent of the putative shell based on the radio continuum filaments and 
our Hi observations is 4°. 4 x 3.°. 4, which converts to 138 pc x 107 pc at 1.8 kpc. If we 
adopt the velocity centers of the highest-velocity clumps given in § 4 (+35 km s" 1 ) as 
an endcap velocity of the shell, its expansion velocity is 55 km s _1 . To derive the total 
Hi mass and kinetic energy of the FVW172.8+1.5, we need to estimate the total mass of 
the expanding shell. This can be done by extrapolating the mass distribution we derive 
for high velocities. Fig. 11 shows the distribution of Hi mass in each 2.2 km s _1 velocity 
interval between t>LSR = +14 and +46 km s" 1 . We assume the receding and approaching 
hemispheres of the expanding shell to be symmetric and centered at vlsr = — 20 km s _1 
. Note that the approaching portion of the shell is not visible because of confusion by 
the Galactic background/foreground Hi emission. While this may not be true, the resulting 
kinetic energy would be roughly correct should the energy injection be spherically symmetric. 
A Gaussian shape is of ten used for the extrapolat i on when fitting the m ass distribution of 



observed Hi shells (e.g. iGiovanelli fc Hayneslll979l ; iKoo et al.lll990l 120061 ). If we adopt this 



Gaussian extrapolation, (the solid line in Fig. 11), we obtain a value of 1.3 x 10 4 ^ 8kpc M Q 
. However, this may overestimate the total mass because for an ideal expanding thin shell 
of uniform density, the mass per unit radial velocity interval is constant. The mass profile 
at the highest velocities could be decreasing due to the velocity dispersion, or the velocity 
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structure, within the shell. However, the mass profile at the central velocities would remain 
constant unless the turbulent velocity in the shell were large enough to be comparable to its 
expansion velocity. 

Here, we develop a simple, but physical, shell model that yields a flat mass profile for the 
central velocities. We consider a spherical shell with a radius R =2°, a thickness AR = 10', 
an expansion velocity v cxp = 55 km s~ x , and a dispersion velocity v a = 5.5 km s _1 . These 
parameters are derived from the observed Hi shell. We assume that the gas density within 
the shell is constant. For the radial profile of expansion velocity within the shell, we consider 
two different cases. The first is that the expansion velocity is constant with radius. The 
second is that the expansion velocity varies linearly with radius. Linearly decreasing or 
increasing velocity stru ctures within the shell are derived by th eoretical studies for different 
ages of radiative shells ( ISlavin fc Coxlll992l ; iBlondin et al.lll998l ). Here, we choose a velocity 
structure for which the expansion velocity at the outer radius drops to 50% of that at the 
inner radius. 

To fit the derived mass distribution, we adopt a three dimensional cube composed of 
451 3 pixels. The z axis is along the line of sight (LOS). The number density of Hi atoms for 
a pixel centered at r = (x 2 + y 2 + z 2 ) 1 ^ 2 is n(r), which is assumed to be constant in our model. 
This pixel has an LOS velocity vi os (x, y, z) = v exp (r)z/r where v exp (r) is the radial expansion 
velocity at r. The Hi 21-cm emission line from this pixel has a Gaussian distribution with 
the LOS velocity dispersion of v a (=5.5 km s _1 ). Then the column density, N(x,y,v), at 
(x, y) over the LSR velocity interval from v to v + Av and the Hi mass, M(t>), between v 
and v + Av, are given as follows: 



N(x,y,v) 



n 



TTVt 




v+Av 



exp 



~{v' - vios(x,y,z) f 
2v 2 



dv'dz, 



M(v) 




N(x, y, v)dxdy. 



The dotted line in Fig.[TT]is the best fit for the mass distribution of a shell with a constant 
expansion velocity profile. Here, the only fitting parameter is the density of the shell. The 
velocity range, where the mass of the observed Hi shell smoothly decreases, is wider than 
the velocity range due to velocity dispersion, so the observed Hi mass distribution cannot 
be explained by velocity dispersion alone. The dashed line shows the best fit for the mass 
distribution of a shell with an expansion velocity profile linearly decreasing towards larger 
radii. This dashed line fits the observed mass better than the dotted line. Variation of the 
expansion velocity with radius dilutes the rectangular mass distribution at high velocities, 
but still yields a constant mass distribution at lower velocities, from vlsr ~ ~ 50 to +10 km 
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s -1 . Note that the basic shape of the mass profile is determined by the difference between the 
minimum and maximum expansion velocities. The detailed velocity and density structures 
only affect the shape of the decreasing portion of the mass profile. We will study the 
mass profiles resulting from various velocity and density structures in a forthcoming paper. 
The estimated total Hi mass of the shell adopting a linearly decreasing velocity profile is 
5900g^ 8kpc M . We will use this value as the characteristic mass of the shell. Note that the 
total mass derived from the Gaussian fit is about twice this. The kinetic energy of the shell, 
taking into account the He abundance of 10% by number, is 0.25 x lO 51 ^ 8kpc ergs. The 
derived physical parameters of the model Hi / continuum shell are presented in Table [5j 



5.1.2. The Origin of FVW172.8+1.5 and the X-rays 



We can think of two possibilites for the origin of FVW172.8+1.5; stellar winds from OB 
stars or SN explosions. HII regions can also produce expanding Hi shells, but their expansion 
velocities and energies are much smaller than those of FVW 172.8+1.5. Furthermore, the 
Hll regions are mostly located at the boundary of the structure, so that they do not appear 
to be major sources of the 3°-4°-sized shell. We first consider whether the stellar winds from 
OB stars can produce FVW172. 8+ 1.5. There ar e more than a dozen O stars in this area 
(IGarmany et al.lll982l ; iMaiz-Apellaniz et al.M2004i ) with eight of them projected against the 
shell (Figs. 5 and 9). We summarize the parameters of these stars in Table 6. Their spectral 
types range from 09.5 to 07, and four are associated with the Sharpless HII regions. The 
other four coul d be possible candidates for the origin of the Hi/continuum shell. If we use 
the relation of iHowarth fc Prinjal (119891 ) and a wind speed of ~ 2,000 km s _1 , then the 
typical wind luminosities of 07-09 main-sequence stars are (5-0.8) xlO 35 erg s -1 . For giant 
O stars, this would be larger by a factor of two. However, the required wind luminosity to 
expla in the parameters of the shell is L w = 77Ek/ (9R/v exp ) = 6 x 10 37 erg s -1 ({Weaver et al. 



19771 ). which is two orders of magnitude greater than the available wind luminosities. Even 
if there are groups of stars in this area (see next section), their contribution to the wind 
luminosity will not be significant because they are no brighter than 09. 

We thus consider that FVW172.8+1.5 is likely to be the result of an SN explosion. 
Certainly, the derived parameters are fully consistent with this interpretation. The required 



SN explosion energy is -Esn — 6.8 x 10 43 nQ 16 .R 



3.16„ 



y exp 



1.35 v 0.161 

A, 



ergs, where Uq is the ambient 



H density in cm 3 , R s is the radius of the shell in pc, t> eX p is the expansion velocity in km 



s , and x is the metalicity in units of the solar metalicity ( ICioffi et al.l ll988). For the above 



parameters, -Esn = 1.3 x 10 51 ergs, where the solar metalicity has been adopted. This value 



of Eg jy is close to the canonical outburst energy of a single SN explosion of 1 x 10 51 ergs 
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The age of the shell, assuming this to be in the radiative stage of SNR evolution, is t — 
0.3 R s /v exp = 0.33 Myr. The SNR interpretation is also consistent with the presence of 
X-rays in this region. The thermal energy of the X-ray emitting gas 3 x 10 50 ergs, which is 
25% of the initial explosion e nergy and close to that expected for a 0.33 Myr-old radiative 
SNR (e.g., ICioffi et al.l Il988h . The apparent segregation between the fast-moving Hi gas 
and the X-ray emitting gas could be due to a non- uniform and/or inhomogeneous density 
distribution in the ambient medium. It is also possible that more than one SN exploded 
within the past few 10 5 yr. We have checked for known pulsars in this area, but none are 
known within the boundary of the putative shell. The nearest pulsar, PSR J0540+32 at 
{£,b) = (176?80, +0?72), is ~ 4° distant from the center of the Hi/continuum feature, and 
lies at a distance of ~6.2 kpc. PSR J0540+32 is thus unlikely to be associated with the 
progenitor of FVW172.8+1.5. 



5.2. Speculations on the Star-Formation History of G173+1.5 



The HII complex G173+1.5 is composed of five Sharpless Hll regions S231-235. Of 
these, all but S234 are clustered around a giant molecular cloud (Fig. 9). There are at least 
fourteen embedded star-forming clusters that are 3-5 Myr old around these four HII regions 
(IKirsanova et al.l 120081 ; iDewangan &: Anandaraol 1201 It ICamargo et al.l l201ll . and references 
therein). S234, located 3° below the other HII re gions, also has an associated cluster with 
young stellar objects having an age of ~ 3 Myr (I Jose et al.l 120081 ) . It has been proposed 
that the formation of some of these star-forming clusters could have been triggered by the 
first-generation stars in these HII regions. 

The presence of an old SNR filling the HII complex G173+1.5 has an interesting impli- 
cation. The five Sharpless HII regions are located within the boundary of the Hi/continuum 
structure, suggesting that their formation could have been triggered by SN explosions, stellar 



winds, or expanding HII regions from a previous generation of stars (See lElmegreenl 11998 



for a review of triggered star formation). The age of the SNR for which we have found 
evidence is, however, only 0.33 Myr, so that it is not possible that the current expanding 
shell triggered the formation of these HII regions. Instead it could be the cluster to which 
the SN progenitor belonged that triggered the formation of these HII regions. There are 
two stellar associations in this area; Aur OBI at a distance of 1.3 kpc and Aur OB2 at 3.2 
kpc ([Humphreys! Il978l ). We suspect that Aur OBI, which is scattered over a wide area, 
(£ = 168.°1-178.°1, b = — 7.°4-4.°2), could be composed of several associations at slightly 
different distances. The dist ance moduli of Au r OBI stars indeed range from 9.69 to 11.67 
or 0.9 to 2.2 kpc in distance jHumphrevslll978h . HD 36483 (09IV) and HD 35921 (09.5111) 
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that are relatively close to the central area of the Hi/continuum structure also belong to Aur 
OBI. If one, or both, of these also belonged to the appropriate cluster, then since the main- 
sequence lifetime of a 20 M Q 09 star is ~ 8 Myr (ISchaller et al.lll992l ). the star formation 
in this area began about £ 8 Myr ago. In this scenario, the first generation stars triggered 
the formation of a second generation that are currently exciting the Sharpless HII regions 
(S231-S235), and we are now observing the formation of third generation stars around these 
HII regions. The HII complex G173+1.5 could be another good example of sequential star 
formation over several stellar associations. 



6. Summary 

We have made high resolution Hi 21-cm line observations of a 7?5 x 4?5 area centered 
at (£,b) = (173°, —1.5°) using the ALFA 7-beam receiver mounted on the Arecibo 305 m 
telescope. The data cube provides a detailed view of the fast-moving Hi gas at velocities 
forbidden by Galactic rotation in this area. FVW 172.8+1.5 is a feature that appears as a 
faint wing in low-resolution surveys. By comparison with radio continuum, X-ray, and Ha 
maps, we find evidence that FVW 172.8+1.5 is likely to be an old SNR within a large HII 
complex situated in this area. Our results are summarized as follows: 

(1) The Hi emission associated with FVW172.8+1.5 is visible from wlsr =+19.5 to +50.4 km 
s _1 . At the highest velocities, the Hi appears to be separated into two concentrations. At 
lower velocities a thin, filamentary feature appears surrounding these two concentrations. 
The increase of the surface area and the amount of Hi gas with decreasing velocity suggests 
that it is an expanding shell. We attribute the complex morphology to the non-uniformity 
and inhomogeneity of the ambient medium. 

(2) The high-velocity Hi emission features are confined inside the radio continuum filaments 
associated with the large HII complex G173+1.5 and show a very good spatial correlation 
with them. This strongly suggests that it has an origin within the complex. Adopting 
a distance of 1.8 kpc, and the systemic velocity of —20 km s _1 of the HII complex to be 
appropriate for the Hi shell, we derive physical parameters for the shell. These yield a kinetic 
energy of 2.5 x 10 50 ergs. This large kinetic energy and the fast-expansion velocity implies 
a SN origin. 

(3) We investigated the nature of the filaments in HII complex G173+1.5 using multi- 
wavelength data. Our analysis of the radio continuum spectrum suggests that the radio 
filaments have both thermal and non-thermal components. A detailed comparison with the 
associated Ha filament supports this conclusion. We also find X-ray emitting hot gas inside 
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the complex with an estimated thermal energy of 3 x 10 50 ergs. This associated hot gas of 
comparable thermal energy to the kinetic energy of expansion, and the probable presence of 
non-thermal radio filaments, supports the conclusion that the HV Hi gas is of SN origin. 

(4) From our analysis, we conclude that FVW 172.8+1.5 is likely to be an old (~ 0.33 Myr) 
SNR produced inside the HII complex G173+1.5. We propose that the stellar association to 
which the progenitor belonged could have triggered the formation of the OB stars currently 
exciting the HII regions. The HII complex G173+1.5 appears to be a potential example of 
sequential star formation over several stellar associations. 
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Table 1. Sharpless Hll regions in the observed field 





r a t) pc n?nnn N i 




^LSR 


^LSR 


D i a n pp 1 ^ 


± 'i Clll-LC 


( h m M 


(') 


( V m ~ ^ ^ 

I -tVlll o / 


firm t; \ 




S229 


5 16.3, +34 27 


65 


+4.4 


+6.7 


0.51 


S230 


5 22.5, +34 08 


300 


+0.0 






S231 


5 39.3, +35 56 


12 


-17.5 


-18.1 


2.3 


S232 


5 42.5, +36 12 


40 


-13.7 


-23.0 


1.0 


S233 


5 38.7, +35 48 


2 


-14.5 


-18.4 




S234 


5 28.1, +34 26 


12 


-14.3 


-13.4 


2.3 


S235 


5 41.0, +35 51 


10 


-25.7 


-18.8 


1.6 


S236 


5 22.6, +33 22 


55 


-3.8 


-7.2 


3.2 


S237 


5 31.4, +34 17 


7 


+1.4 


-4.3 


1.8 



Optical size (|Sharplesslll959l ) 



b Ha emission line radial velocity ( Fich et alJll990i ) 



c CO radial velocity (jBlitz et al.lll982l ). For S229, the association is uncertain. 

d Distance dete rmined by spectro photometry of the central exciting star. See 
the references in lBlitz et al.l (|1982l ). 



Table 2. Parameters of the continuum surveys used for spectral analysis 





Frequency 






RMS 




Name 


(MHz) 


Resolution 




(mK) 


References 


WENSS 


325 


54 "x54 "esc 


8 


2300 


Reneelink et al. ( 1997) 


CGPS 


408 


2 .'8 x 2'. 8 esc 


S 


950 


Tavlor et al. (2003) 


CGPS 


1420 


49 "x49 "esc 


S 


68 


Tavlor et al. C2003) 


Effelsberg 11-cm 


2695 


4'.3 




20 


Fiirst et al. (1990) 
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Table 3. Integrated radio continuum flux densities of Filament A and S232 





325 MHz 


408 MHz 


1420 MHz 


2695 MHz 


Spectral index 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


(S oc u~ a ) 


Filament A 


1.5 ±0.17 


0.84 ±0.14 


0.87 ±0.03 


0.73 ± 0.05 


±0.23 ± 0.06 


S232 


1.8 ±0.07 


1.9 ±0.10 


1.9 ±0.10 


2.2 ±0.05 


-0.09 ± 0.02 



Note. — The areas used for the flux density measurements of the filament A and 
S232 are indicated in Fig. [71 





Table 4. 


Spectral indices from T-T plots 






325-1420 MHz 


325-2695 MHz 408-1420 MHz 


408-2695 MHz 


Filament A 


±0.50 ±0.11 


±0.32 ±0.06 ±0.07 ±0.10 


-0.04 ±0.08 


S232 


±0.09 ± 0.07 


-0.00 ± 0.04 ±0.06 ± 0.04 


-0.04 ± 0.02 



Table 5. Parameters for an expanding shell 



Parameter Estimated value 



Center (a, 5) 


05 h 33 m 29 s ,±35°50 / 52 // 


(£,b) 


172°.78,±1°.51 




Mean radius 


69 x 54di. 8 kpc pc 




Expansion velocity 


55 ± 5 km s _1 




Total Hi mass 


5900 ± 770 4 8kpc 


M 


Kinetic energy 


0.25 ± 0.05 x 10 51 


4 8 kpc er g S 


Initial ISM density 


0.25 ±0.03 d^ 8kpc 


cm -3 


Age 


0.33 ± 0.03 Myr 





Note. - The errors are la errors. For the esti- 
mation of the initial ISM density and age, we use the 
geometrical average of the mean radius (61 pc). 
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Table 6. O-type stars within the HII complex G173+1.5 



RA, Dec (J2000) Spectral Distance Associated 
Name ( h m, ° ' ) Type ( kpc ) Hn Region 



HD 36483 


5 


33.8, 


+36 


28 


09IV 


1.3 




HD 35921 


5 


29.8, 


+35 


20 


09.5111 


1.3 




LS V+35 24 


5 


39.7, 


+35 


53 


09V 


2.1 


S231 


BD +35 1201 


5 


41.0, 


+35 


50 


09.5V 


1.6 


S235 


HD 37737 


5 


42.4, 


+36 


09 


09.5111 


1.3 


S232 


BD +34 1058 


5 


28.4, 


+34 


40 


08 


3.2 




HD 35619 


5 


27.5, 


+34 


47 


07V 


3.2 




BD +34 1054 


5 


28.1, 


+34 


25 


09.5V 


2.0 


S234 



Note. — The O stars ([Garmanv et al.lll982l ; lMaiz-Apellaniz et al.ll2004l ) 
located within a radius of 2°5 from the center of the Hi/continuum struc- 
ture are listed in order of angular distance from the center. See also Fig. [5j 
The distance of a star is computed from its distance modulus, ass uming 
an interstellar absorption of A(V)=3E(B-V) (jGarmany et al.lll982l ). 
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RA (J2000) RA (J2000) 

Fig. 1. — Hi channel images of the observed field. Each image is integrated over a velocity 
interval of 5.1 km s _1 . The central LSR velocity is given at the bottom left of each image. 
The numbers in the brackets are the minimum and maximum brightness temperatures (K) 
of the grey scale which is linear. The black cross indicates the approximate center of the 
Hi structure. The dotted ellipse on the +24.7 km s -1 frame marks the boundary of the 
proposed shell as adopted in this paper. 
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5 h 55 m 50 45 40 35 30 25 20 

RA (J2000) 



Fig. 2. — A three-color image showing the Arecibo Hi emission from FVW172.8+1.5. Red, 
green, and blue represent the images integrated over LSR velocities of +45 to +35, +35 
to +25, and +25 to +20 km s™ 1 respectively. Effelsberg 11-cm radio continuum contours 
(jFiirst et al.lll990l ) are overplotted. Contour levels are 30, 100, 200 mK in brightness tem- 
perature. 
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Fig. 3. — The Hi peak brightness temperature over the range of Galactic longitude, 171. °8 
to 173. °7, as a function of Galactic latitude. In order to show Hi features at lower velocities, 
the intensity scale for brightness temperatures in excess of 5 K is greatly expanded. The 
grey-scale ranges are shown by the bar in the upper right of the figure; values for brightness 
temperatures < 5 K are marked below the bar, with those for > 5 K are marked above. The 
merit of a peak intensity image is that it enhances bright structures that would not be seen 
in an averaged image. 




Fig. 4. — Sample line profiles in the area around (5 h 34. 5 m , +36° 03'). The relative positions 
are shown in each frame. Dashed and thick solid lines show the original and background- 
subtracted line profiles, respectively. The thin solid lines show the best fit of Gaussian profiles 
to the background-subtracted profiles. The central velocities and velocity dispersions of the 
Gaussians in km s _1 are given in each frame. 
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5 h 50 m 40 30 20 10 

RA (J2000) 



Fig. 5. — A grey-scale plot of the Effelsberg 11-cm radio continuum in the direction of 
FVW172.8+1.5. Sharpless Hll regions are marked by blue circles, while Galactic O stars are 
marked by red star symbols. The region observed at Arecibo (Fig. 1) is shown by green lines. 
The black cross and the dotted ellipse are the approximate center and boundary adopted 
for the shell, respectively. The continuum filaments defining the Hll complex are marked by 
arrows. The scale bar in the bottom left represents 50 pc at 1.8 kpc. 
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RA (J2000) 

Fig. 6. — An Hi map integrated from vlsr — — 28 to —25 km s" 1 , with the 11-cm continuum 
contours overlaid. An HI filament along the continuum filament A is visible (cf. Fig. 5), 
although the correlation is not perfect. 
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174 173 172 

GLON (deg) 



Fig. 7.— WENSS 325 MHz (top left), CGPS 408 MHz (top right), CGPS 1420 MHz (bottom 
left), and Effelsberg 2695 MHz (bottom right) radio continuum images of the northern part 
of the FVW173+1.5 field. The areas used to derive the total flux density and to plot T- 
T diagrams of radio continuum filament A and the Hn region, S232, are marked by thick 
dotted and solid boxes, respectively. The neighboring areas used for the subtraction of the 
large scale background are marked by thin dotted and solid boxes. The white holes in the 
images are where point sources have been subtracted. 
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Fig. 8. — (Left) Flux densities of the Hll region S232 (open circles) and the continuum 
filament A (filled circles). The fitted spectra, with a (S oc z/~ Q ) = -0.09 ±0.02 and +0.23 + 
0.06 for the Hll region and filament respectively, are overplotted. (Right) The T-T plot for 
the area of filament A using the WENSS 325 MHz and Effelsberg 2695 MHz data. The best 
linear fit with a = +0.32 ± 0.06 is overplotted. 
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Fig. 9. — (Top left) The ROSAT R45 (3/4 keV) band diffuse X-ray contour map overplotted 
on the Effelsberg 11-cm image. The contour levels are 70, 85, and 100 in counts s _1 arcmin~ 2 . 
(Top right) As for the top left frame, but for the ROSAT R67 (1.5 keV) band. The contour 
levels are 100, 120, 140, and 160 in counts s -1 arcmin -2 . (Bottom left) The Ha map (6' 
FWHM resolution) composite from t he Virginia Tech Spectral line Survey (VTSS) and the 
Wisconsin H-Alpha Mapper survey ( IFinkbeinerl 120031 ) . The scale for the image is 40 - 70 
Rayleigh, but a histogram equalization method is performed to enhance the contrast of the 
H-alpha filaments. (Bottom right) CO J=l-0 intensity image integrated over wlsr — — 25 to 
— 15 km s _1 . The scale for the image is - 0.3 K, and is linear. The blue ellipses in solid 
and dotted lines overplotted on the top two images are the area used to measure the X-ray 
counts of the shell and the backgrounds, respectively. The O stars are marked by red stars. 
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Fig. 10. — (Left) An enlarged view of the filament A in the VTSS Ha image with the 11-cm 
radio continuum contours superposed. (Right) One-dimensional intensity distributions of 
Ha (dotted lines) and 11-cm continuum (solid lines) emission along Galactic latitudes at 
several Galactic longitudes. Ha and continuum intensities are normalized by 20 R and 0.1 
K, respectively. 
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Fig. 11. — The mass distribution for the putative Hi shell assuming a distance of 1.8 kpc. 
The filled circles are the derived Hi masses for each 2.2 km s _1 velocity interval. The error 
bars show the standard deviation for the background area. The solid line shows a Gaussian 
fit to the data, while the dotted line shows the best fit for the mass distribution of a shell 
with a constant v cxp =55 km s _1 and v a = 5.5 km s _1 , and the dashed line shows the best 
fit for the mass distribution of a shell with a velocity dispersion v a = 5.5 km s" 1 , but whose 
expansion velocity at the outer radius of the shell drops linearly to 50% of the expansion 
velocity at the inner radius. 



